clc;
clear;

gauss_points = [(1-sqrt(3/5))/2, (1+sqrt(3/5))*(1-sqrt(3/5))/4;
    (1-sqrt(3/5))/2, (1+sqrt(3/5))/4;
    (1-sqrt(3/5))/2, (1+sqrt(3/5))*(1+sqrt(3/5))/4;
    1/2, (1-sqrt(3/5))/4;
    1/2, 1/4;
    1/2, (1+sqrt(3/5))/4;
    (1+sqrt(3/5))/2, (1-sqrt(3/5))*(1-sqrt(3/5))/4;
    (1+sqrt(3/5))/2, (1-sqrt(3/5))/4;
    (1+sqrt(3/5))/2, (1-sqrt(3/5))*(1+sqrt(3/5))/4];
gauss_weights = [100/324*(1+sqrt(3/5))/8, 5/81*(1+sqrt(3/5)), 100/324*(1+sqrt(3/5))/8, 5/81, 8/81, 5/81, 100/324*(1-sqrt(3/5))/8, 5/81*(1-sqrt(3/5)), 100/324*(1-sqrt(3/5))/8];

vpa(gauss_points)

